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Abstract. 

A new representation, which does not contain the third-order derivatives of 
the coordinates, of the exact Mathisson-Papapetrou-Dixon equations, describing the 
motion of a spinning test particle, is obtained under the assumption of the Mathisson- 
Pirani condition in a Kerr background. For this purpose the integrals of energy 
and angular momentum of the spinning particle as well as a differential relationship 
following from the Mathisson-Papapetrou-Dixon equations are used. The form of these 
equations is adapted for their computer integration with the aim to investigate the 
influence of the spin-curvature interaction on the particle's behavior in the gravitational 
field without restrictions on its velocity and spin orientation. Some numerical examples 
for a Schwarzschild background are presented. 



PACS numbers: 04.20.-q, 95.30. Sf 
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1. Introduction 

In general relativity two main approaches have been developed for the description of 
spinning particle behavior in a gravitational field. Chronologically, the first one was 
initiated in 1929 when the usual Dirac equation was generalized for curved space-time 
[1]. The second, the pure classical (non-quantum) approach, was proposed in 1937 [2]. 
Later it was shown that in a certain sense the equations from [2] follow from the general 
relativistic Dirac equation as a classical approximation [3]. 

The focus of this paper is on the equations of motion of the classical spinning 
particle, which after [2] were obtained in [4] and in many later papers by different 
methods. These equations can be written as 

DS"" „ DS'''' , DS"" 

+ u"u^— u'^u^—— = 0, 2 

ds ds ds 

where = dx^ /ds is the particle's 4- velocity, S"'' is the tensor of spin, m and D/ds 
are, respectively, the mass and the covariant derivative with respect to the particle's 
proper time s; R\p^ is the Riemann curvature tensor (units c — G — 1 are used); here 
and in the following, latin indices run 1, 2, 3 and greek indices 1, 2, 3, 4; the signature 
of the metric (-,-,-,-!-) is chosen. 

Equations (1), (2) were generalized in [5] for the higher multipolcs of the test 
particles and now the set (1), (2) is known as the Mathisson-Papapctrou-Dixon (MPD) 
equations. Multipolar equations of motion for extended test bodies in general relativity 
were considered in a recent paper [6] and in this context the importance of the seminal 
work [2] was pointed out. 

The first effects of the spin-gravity interaction following from (l)-(2) were 
considered in [7] for the Schwarzschild field. According to [7] and in many further 
publications (a list of them is presented, for example, in [8, 9]) the infiuence of spin on 
the particle's trajectory is negligible small for practical registrations. However, in this 
sense much more realistic are the effects connected with spin precession [10]. 

An interesting point has been elucidated in [11] concerning the possibility of a 
static position of a spinning particle outside the equatorial plane of the Kerr source 
of the gravitational field, on its axis of rotation. In spite of the conclusion that such 
a situation is not allowed by the MPD equations [11], this question stimulated the 
investigations of possibilities of some non-static (dynamical) effects connected with the 
particle's motion relative to a Schwarzschild or Kerr mass outside the equatorial plane 
[12]. Then it was shown that spinning particles moving with relativistic velocity can 
significantly deviate from geodesies [13, 14]. 

Some papers are devoted to the investigation of equilibrium conditions of spinning 
test particles in the Schwarzschild-de Sitter and Kerr-de Sitter space-times [15], as a 
development of studying the stationary equilibrium positions of charged particles in 
Reissner-Nordstrom and Kerr-Newman space-times [16]. 
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While investigating the solutions of equations (1), (2), it is necessary to add a 
supplementary condition in order to choose an appropriate trajectory of the particle's 
center of mass. Most often conditions [2, 17] 

S^^'u^ = (3) 

or [5, 18] 

S^'^P^ = (4) 
are used, where 

P'^ = mu" + ux^^ (5) 

ds 

is the 4-momentum. The condition for a spinning test particle 

^ ^ £ « 1 (6) 
mr 

must be taken into account as well [11], where l-S'oj = const is the absolute value of spin, 
r is the characteristic length scale of the background space-time (in particular, for the 
Kerr metric r is the radial coordinate), and Sq is determined by the relationship 

-^0 = Is.^S'^''. (7) 

In general, the solutions of equations (1), (2) under conditions (3) and (4) are 
different. However, in the post-Newtonian approximation these solutions coincide with 
high accuracy, just as in some other cases [19, 20]. Therefore, instead of exact MPD 
equations (1) their linear spin approximation 

D A 1 
ds 2 

is often considered. In this approximation condition (4) coincides with (3) (by condition 
(3) m in equations (1) is a constant quantity). 

According to [21] for a massless spinning particle, which moves with the velocity of 
light, the appropriate condition is (3). The question is which condition is adequate for 
the motions of a spinning particle with the nonzero mass if its velocity is close to the 
velocity of light? To answer this question it is necessary to analyze the corresponding 
solutions of the exact MPD equations (1), (2) both at condition (3) and (4). 

The main purpose of this paper is to consider the exact MPD equations under 
condition (3) in a Kerr metric. Due to the symmetry of this metric equations (1), 
(2) have the constants of motion: the particle's energy E and the projection of its 
angular momentum [22-25]. It is known that in the case of the geodesic equations 
the analogous constants of motion were effectively used for analyzing possible orbits of 
a spinless particle in a Kerr space-time [26, 27]. Namely, by the constants of energy 
and angular momentum the standard form of the geodesic equations, which are the 
differential equations of the second-order by the coordinates, can be reduced to the 
differential equations of the first order. Naturally, it is interesting to apply a similar 
procedure to the exact MPD equations. However, in contrast to the geodesic equations. 



- -^u^S^^R\^, (8) 
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the exact MPD equations at the condition (3) contain the third derivatives of the 
coordinates [28, 29]. Therefore, the apphcation of this procedure to the exact MPD 
equations is significant. 

In this paper in order to obtain a full set of the MPD equations at condition (3), 
without the third derivatives of the coordinates, some differential relationship following 
from equations (1), (2) is used. We present this relationship in section 2 in general form, 
for any metric. Its concrete form in the Kerr metric, together with the expressions for 
E and is used in section 3 and the full set of the differential equations for the 
dimensionless quantities connected with the particle's coordinates, velocity and spin is 
described. The explicit form of these equations are written in the Appendix. We analyze 
the relationship between u'^ and at condition (4) in section 4. Section 5 is devoted 
to some numerical examples. We conclude in section 6. 

It is important to note that the very condition (3) arose in a natural fashion in the 
course of its derivation by different methods [30-32]. Therefore, it is of importance to 
obtain a representation of the exact MPD equations at this condition in the Kerr metric 
convenient for their further computer integration. 

We point out that the integrals of energy and angular momentum of the MPD 
equations in a Kerr space-time were effectively used for different purposes in [8, 22-25, 
33-38] at condition (4). 



2. A relationship following from equations (l)-(3) 

In addition to the antisymmetric tensor S'^" in many papers the 4- vector of spin sx is 
used as well, where by definition 

sx = l^^ex^.^^S'"' (9) 

and g is the determinant of the metric tensor, ex^iy^j is the Levi-Civita symbol. It 
follows from (7), (9) that sxs^ = Sq and at condition (3) we have Sxu^ = (other 
useful relationships with sx following from MPD equations at different supplementary 
condition can be found, for example, in [9]). 

The set of equations (2) contains three independent differential equations and in 
(3) we have three independent algebraic relationships between S^'^ and u^. By (3) the 
components 5"*^ can be expressed through S'^'^: 

So, using (10) the components 5"^ can be eliminated both from equations (2) and (1). 
That is, in further consideration one can "forgot" about supplementary condition (3) 
and deal with the three independent components 5'*'^. However, it is appear that more 
convenient form of equations (1), (2) is not for S"^^ but for another 3-component value 
Si which is connected with S^^ by the simple relationship 
1 
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where Eiki is the spatial Levi-Civita symboL For example, it is not difficult to check that 
three independent equations of set (2) in terms of Si can be written as 

U^Si + 2(u[4Ui] - u''Upr'^^^Ui])Sku'' + 2SnK[4Ui]U'' = 0, (12) 

where a dot denotes usual differentiation with respect to the proper time s, and square 
brackets denote antisymmetrization of indices; T^^ are the Christoffel symbols. 

The simple calculation shows that the 3-component value Si has the 3-vector 
properties relative to the coordinate transformations of the partial form — 
x^ {x^ , x"^ , x^) , x'^ = x"^ and in this special sense Si can be called as a 3-vector. (By 
the way, in the context of equations (1), (2) firstly the 3-vector of spin was used in [7] 
with the notation S = (S'^'^, S'^\ S^'^)). By equations (9)-(ll) the relationship between 
Si and s\ is 

Si = -Si + —S4. (13) 
U4 

Let us consider the first three equations of the subset (1) with the indexes A = 1, 2, 3. 
Multiplying these equations by Si, S2, S^ correspondingly and taking into account 
(9)-(ll) we get 

mS,^ = --u-S^'^S.R^^^ (14) 

(here the covariant derivative Du^/ds remains 4D, i.e., is determined according to the 
Christoffel connection of the 4-dimensional space-time). We stress that in contrast to 
the each equation from set (1), which contain the third derivatives of the coordinates, 
relationship (14) does not have these derivatives. 

Relationship (14) is an analog of relationship (21) from [9] where the spin 4- vector 
s\ is used. 



3. On set of exact MPD equations with constants of motion E, for the 
Kerr metric 



In the Boyer-Lindquist coordinates x^ = r, x^ = 6, x^ = ip, x'^ = t the non-zero 
components of the Kerr metric tensor are 

9u = 922 = -p , 



where 



2Mra'^ 

7 



«?33 = - ( + + — sin^ 9 ] sin^ 9, 



2Mra .3^ ^ 2Mr ,^ , 

534 = — 5— sm 9, 5^44 = 1 5-, (15) 



p2 ^ + cos^ ^, A = r^-2Mr + a^ < ^ < tt. 
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In this coordinates the constant values of the particle's energy E and the projection 
of its angular momentum can be written as [22-25] 

E^P,-^g,,,,S^\ (16) 

J, = -P3 + ^^3M,.^'^^ (17) 

It is convenient to use the dimensionless quantities t/i connected with the particle's 

coordinates by definition 

r t 
Vi ^j^, y2 = 0, Vz = y4 = (18) 

as well as the quantities connected with its 4- velocity 

Vb = u\ 176 = Vi = Mu\ ys = (19) 

and the spin components [14] 

•^i S2 S3 , , 

(We underline that in the following, in sections 4, 5 and in the Appendix, the notation 
means that the quantity number i from the set of eleven quantities (18)-(20) is to 
the k-th power.) In addition, we introduce the dimensionless quantities connected with 
the particle's proper time s and the constants of motion E , J/. 

^ = 4-.^ E = -, J = A7- (21) 
M m mM ^ ' 

Quantities (18), (19) satisfy the four simple equations 

yi = ys, 2/2 = ye, ys = y?, y4 = ys, (22) 

here and in the following a dot denotes the usual derivative with respect to x. 

Now we point out other seven nontrivial first-order differential equations for the 11 
functions y^. Namely, the first of them follows directly from equation (14). The second 
is a result of the covariant differentiation of the normalization condition Uyvy — 1, that 
is 

= 0. (23) 

ds 

The third and fourth equations follow from (16) and (17) correspondingly if condition 
(3) is taken into account. Finally, the last three equations for y^ follow directly from (12). 
This set of the seven equations is presented in the Appendix. Equations (A.3)-(A.9) 
together with the four equations from (22) are the full set of the exact MPD equations 
which describe most general motions of a spinning particle in the Kerr gravitational 
field without any restrictions on its velocity and spin orientation. We stress that the 
two equations following from (16) and (17) (see (A. 5), (A. 6)) contain the quantities E 
and J as the parameters proportional to the particle's energy and angular momentum 
according to notation (21). 
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Now we recall some features of the solutions of the exact MPD equations under 
condition (3). It is known that in the Minkowski space-time equations (l)-(3) have, 
in addition to usual solutions describing the straight worldlines, a set of solutions 
describing the oscillatory (helical) worldlines [28, 29]. The physical interpretation of 
these superfluous solutions was proposed by C. MoUer [30]. He pointed out that in 
relativity the position of the center of mass of a rotating body depends on the frame of 
reference, and condition (3) is common for the so-called proper and non-proper centers 
of mass [39] . The usual solution describe the motion of the proper center of mass and 
the hehcal solutions describe the motions of the set of the non-proper centers of mass. 
Naturally, in general relativity, when the gravitational field is present, the exact MPD 
equations (l)-(3) have some superfluous solutions as well. Just to avoid these solutions, 
instead of (3) condition (4) was used in many papers. In contrast to (3), condition (4) 
picks out the unique worldline of a spinning particle in the gravitational field. However, 
the question arises: is this worldline close, in the certain sense, to the usual (non-helical) 
worldhne of equations (1), (2) under condition (3)? It is simple to answer this question 
if the linear spin approximation is valid, because in this case condition (4) practically 
coincides with (3). Whereas another situation cannot be excluded a priori for the high 
particle's velocity. 

Concerning equations (22), (A.3)-(A.9) we stress that by choosing different values of 
E and J for the fixed initial values of yi one can describe the motions of different centers 
of mass. Among the set of the pairs E and J there is the single pair corresponding to the 
proper center of mass. The possible approaches for finding this pair is a separate subject. 
One of them was proposed in [40] where a method of separation of some non-oscillatory 
solutions of the exact MPD equations in the Schwarzschild field was considered. In the 
next section we shall analyze the possibility of using the expressions for E and from 
(16), (17) at condition (4) for the same purpose. Note that under condition (4) the 
concrete values of E and Jz are fully determined by the initial values of the particle 
coordinates and velocity (or momentum), without its acceleration, in contrast to the 
case with condition (3). 

4. Values E and J according to condition (4) 

Let us check the supposition that the single solutions of equations (1), (2) at the 
supplementary condition (4), corresponding to the fixed initial values of the coordinates, 
velocity and spin, is close to those solutions of equations (1), (2) at the condition (3) 
which describe the motion of the proper center of mass with the same initial values. As 
we pointed out in section 1, this assumption is justified for the velocities which are not 
close to the velocity of light. Here we shall consider the situation for any velocity. 

Let us write the main relationships following from the MPD equations at condition 
(4) [22-25]. The mass of a spinning particle m is defined as 
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and m is the integral of motion, that is dm/ds — 0. The quantity is the normahzed 
momentum, where by definition 



m 



(25) 



Sometimes is called the "dynamical 4- velocity", whereas the quantity from (1)- 

(3) is the "kinematical 4- velocity" [8]. As the normalized quantities u'^ and satisfy 
the relationships 

U),v^ = 1, VxV^ = 1. (26) 

There is the important relationship between ti'^ and [22-24]: 

1 



u 



N 



(27) 



where 

A^l + ^Rx^,,S'-SP% (28) 

Now our aim is to consider the explicit form of expression (27) for the concrete 
case of the Schwarzschild metric, for the particle motion in the plane 9 — 7r/2 
when spin is orthogonal to this plane (we use the standard Schwarzschild coordinates 



X — r, 



X 



9, x^ 



ip, x^ = t). Then we have 



u 



0, u^^O, liVO, liVO, 



^12 ^ ^23 ^ ^13 ^ Q 

In addition to (30) by condition (4) we write 



q14 _ -^3 r,13 



(29) 
(30) 

(31) 



Using (7), (29)-(31) and the corresponding expressions for the Riemann tensor in the 
Schwarzschild metric, from (27) we obtain 



u 



u 



u 



3 {V^V - 1) 



0, 



NV^ 1 + 



3M 



C2 



m^A 



(32) 



where M is the mass of a Schwarzschild source. According to (28) we write the expression 
for A as 



A = 1 



(33) 



(the quantity M in (32), (33) is the mass of a Schwarzschild source). Inserting (33) into 
(32) we get 



1 + 



rn^r^ 
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rn^r^ 



1 + 



SlM 
rn^r^ 



As in (6), we note 



mr 



(34) 



(35) 



where according to the condition for a test particle it is necessary £ <^ 1. However, in 
our calculations we shall keep all terms with e. 

The explicit expressions for N we obtain directly from the condition uxu^ = 1 in 
the form 

2 



M 



r I r 



X ( 2-e^—j 



-1/2 



(36) 



Inserting (36) into (34) we obtain the expression for the components through 
(y2^ii2^0): 



,M 



where 



R 



2e' 



.M 



- 2,{u^fe^Mr [2-e 



,M 



-1/2 



(37) 



(38) 



The main feature of relationships (37), (38) is that for the high tangential velocity 
of a spinning particle the values V^,V^, become imaginary. Indeed, if 

1 



\u'\ > 



£^6iWr' ^^^^ 
in (38) we have the square root of the negative value. (As writing (39) we neglect the 
small terms of order e^; all equations in this section before (39) and after (40) are strict 
in e.) Using the notation for the particle's tangential velocity Utang = ^'i*^ by (39) we 
write 



-\/6M' 



(40) 
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According to estimates similar to those which are presented in [14] if r is not much 
greater than M, the velocity vahic of the right-hand side of equation (40) corresponds 
to the particle's highly relativistic Lorcntz 7- factor of order 1/e. 

Probably, this fact that according to (25), (37)-(40) the expressions for the 
components of 4-momentum become imaginary (if m in (24) is real) is an evidence 
that condition (4) cannot be used for the particle's velocity which is very close to the 
velocity of light. However, this point needs some additional consideration. In any case, 
relationships (37)-(40) are of importance for authors which investigate solutions of the 
MPD equations at condition (4). We stress that many papers were devoted to study the 
planar or circular motions of spinning particles in the Schwarzschild or Kerr space-time 
at different supplementary conditions [7-14, 22-25, 33-35, 40, 41 ]. Equations (37)-(40) 
elucidate the new specific features which arise for the highly relativistic motions. 

Another aspect of the connection between the spinning particle's momentum and 
velocity at condition (4) was considered in [22, 23]. It is pointed out in [22] that there 
exist a critical distance of minimum approach of a particle to the Kerr source where 
its velocity becomes space-like with the time-like momentum. Similarly according to 
conclusion from [23] the velocity of the spinning particle become space-like for sufficiently 
large gravitational fields and/or spins. In contrast to [22, 23] expressions (37)-(40) 
describe the situation when the velocity is time-like but the components of momentum 
become imaginary. 

The new result of this section as compare to [22, 23] consists in the conclusion that 
according to (37)-(40) only tangential component of velocity is important in this case, 
not the radial one, although the orbit is not necessarily circular. 

It is interesting to check the possibihty of using the values E and Jz as calculated 
by (37) for computing spinning particle motions by the equations which are described 
in the previous section if the particle's tangential velocity is much less than the critical 
value from the right-hand side of (40) . At condition (4) the constants E and for the 
equatorial motions in the Schwarzschild field can be written as 

E^P, + IquaS'^ = mV, - L,,^,Yls'\ (41) 

Jz = -Pz - ^^33,1-513 = -mV^ - ^^33,1-5^^ (42) 

Using the dimensionless quantities yi as defined in (18), (19), relationship (37) and the 
simple expression for S^^ through 5*0 from (41), (42) we obtain 



E^R 



(43) 
(44) 



yij \ yi J \ vi. 

-] + eyws (1- -) (1- — 

vij \ yij \ yi 

where according to (21) we note E = E/m, J = Jz/mM. Then for the fixed initial 
values of the quantities yi,y7,ys in (43), (44) we have the concrete values of E and 
J which can be used for numerical integration of the exact MPD equations. Some 
examples we shall consider in the next section. 
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5. Numerical examples 

Let us consider some solutions of the exact MPD equations (22), (A.3)-(A.9) for the 
equatorial particle's motion in the Schwarzschild field with the values E, J from (43), 
(44). We are interested in the highly relativistic motions when u^^^^g 3> 1 and, at the 
same time, |xitang| is much less than right-hand side of (40). All figures 1-8 correspond 
to the situation when the initial values of the particle coordinates and velocity are given 
as 

r(0) = 2.5M, ip{0) = 0, u\0) = -1, Mu^(0) = 100 (45) 

(for the equatorial motions 9 — 7r/2 identically, i.e., = 0; u'^{Q) is determined from 
the condition Uj^u'' — 1). The small value So/ Mm is equal to 10~^ for figures 1-3 and 
6 X 10~^ for figures 4-8. According to (35) for the spin 3-vector we have Si — 0, Ss — 
and it follows from (12) that S2 — rSo- For comparison, we present the corresponding 
solutions of the geodesic equations with the same initial values of the coordinates and 
velocity. 

Figures 1 and 4 show the dependence r vs. s, with the difference that in figure 1 the 
graph for a spinning particle practicaUy coincides with the corresponding geodesic graph, 
whereas according to figure 4 the graph for a spinning particle reveals the osciUatory 
features. The similar property takes place for the dependence (p vs. s according to figure 
5 (the simple linear dependence (p vs. s for a spinning particle with So/Mm — 10~^ 
is not presented here). Figures 2, 3, 6, 7 illustrate the clear oscillatory regime for the 
radial and angular particle's velocity. We point out that the amphtude of the oscillation 
increases with the value of spin, whereas the frequency of this oscillation decreases. At 
the same time, according to figure 8 the dependence r vs. the coordinate time t does 
not reveal the oscillatory features. It is not strange because dt/ds oscillates similarly to 
dr/ds and dip/ds (the corresponding graph for dt/ds is not presented here for brevity) 
in such a manner that dr / dt is not oscillatory. 

It is known that non-proper centers of mass of a spinning particle oscillate around 
the proper center of mass [28, 29, 39]. Therefore, one can consider the middle lines of 
the corresponding oscillatory lines in figures 2-7 as such that present the motions of the 
proper center of mass. So, when jufH,,,;! is much less than right-hand side of (45), the 
figures above show that equations (22), (A.3)-(A.9) with the relationships (43), (44) 
can be used to describe the spinning particle motions in the Schwarzschild space-time 
in some approximation. 

Because of the specific oscillatory regime the corresponding graphs in figures 1-8 are 
calculated on the limited proper time interval. To illustrate some longer trajectories we 
shall use the expressions for E and J which follows from relationships (34)-(43) of paper 
[40] . We recall that in [40] an approach for separation of some non-oscillatory solutions of 
the exact MPD equations at condition (3) in the equatorial plane of the Schwarzschild 
space-time was developed (these solutions describe the orbits closer to circular with 
r < 3M). So, according to our notation yi — r/M, E — E/m, J — Jz/mM, and 
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Figure 1. Radial coordinate vs. proper time at So /Mm =10 ^. The line for the 
spinning particle practically coincides with the geodesic line at the same initial values 
(45). 



£o = So/ Mm, now we rewrite the corresponding expressions which are presented in 
equations (34)-(43) of [40] as 



E 



1 



1 - 



J = y/^m 



2 

yi 

2 

1 

yi 



1/4 



1/4 

3 

yi 



-1/4 





3 


-1/2 








1 - 


yi 










-3/2 / 


3 




3 






K 


1 

yi 


+ 


yl 




3 


-3/2 


(- 


9 


+ 


15 










yi 






yi 







(46) 



(47) 



(48) 



Figures 9-11 correspond to the situation when the initial values of the particle 
coordinates and velocity are given as 

r(0) = 2.8M, (^(0) = 0, v}{Q) = -0.5 (49) 

{u^ — identically). In addition we put Sq — 10~^, then according to (46) we have 
Mu? fa —1635. For comparison in figures 9-11 the corresponding graphs for the geodesic 
motion are presented. Figures 9-11 illustrate the significant space separation of the 
trajectories of spinning and spinless particles: for the proper time s 1.3 x 10~^M (it 
corresponds to the coordinate time t ~ 25M) the spinless particle falls on the horizon 
surface after less than 0.5 revolution by the angle (/9 about the Schwarzschild mass, 
whereas the position of the spinning particle by the radial coordinate is close to the 
initial value r = 2.8M for one revolution. 

We point out that figures 9-11 are similar to same figures from our paper [14]. 
However, the essential difference must be stressed: all figures in [14] are presented for 
the linear spin approximation of the MPD equations, whereas here we deal with the 
exact equations. 
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Figure 2. Radial velocity vs. proper time for the spinning particle with So/Mm = 
10~^ (solid line) and for the geodesic motion (dashed line). 
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Figure 3. Angular velocity vs. proper time for the spinning particle with So/Mm 
10~^ (solid line) and for the geodesic motion (dashed line). 
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Figure 4. Radial coordinate vs. proper time for the spinning particle with So/ Mm = 
6 X 10~^ (solid line) and for the geodesic motion (dashed line). 
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Figure 5. Angle tp vs. proper time for the spinning particle with So/Mm = 6 x 10 
(solid line) and for the geodesic motion (dashed line). 
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Figure 6. Radial velocity vs. proper time for the spinning particle with So/Mm = 
6 X 10~^ (solid line) and for the geodesic motion (dashed line). 




0.0000 0.0002 0.0004 0.0006 0.0008 0.001 0.0012 
S (M) 



Figure 7. Angular velocity vs. proper time for the spinning particle with So /Mm = 
6 X 10~^ (solid line) and for the geodesic motion (dashed line). 
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Figure 9. Radial coordinate vs. proper time for the spinning particle with So/Mm = 
10~^ (solid line) and for the geodesic motion with the same initial values of the 
coordinates and velocity (dashed line). 



6. Conclusions 



In this paper we obtained the representation of the exact MPD equations at 
supplementary condition (3) for the Kerr space-time by using the constants of the 
particle's motion, the energy and angular momentum, together with the differential 
consequence of these equations (14). The full set of the corresponding 11 first-order 
differential equations is presented in (22), (A.3)-(A.9). The computer integration of 
these equations is performed for more simple case of the Schwarzschild space-time. The 
possibiUty using expressions (43), (44) to describe motions of a spinning particle in 
this space-time is considered. Naturally, it is not obvious a priori that the constants 
of motion which are calculated assuming condition (4) can be used for the correct 
description of the particle motions by the equations which were derived using condition 
(3). According to figures 1-8 it is possible in some approximation. For a more exact 
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Figure 11. Trajectories in the polar coordinates of the spinning particle with 
So /Mm = 10~^ (solid line) and spinless particle (dashed line) with the same initial 
values of the coordinates and velocity. The circle with the radius 2 corresponds to the 
horizon line. 



description the expressions for E, J from (43), (44) must be corrected. To get these 
improved values of E and J one can use of computer search. The suitable values of the 
pair E, J must give the solutions without large amplitude oscillations. In some cases 
instead of (43), (44) it is possible to use expressions (47), (48) (figures 9-11). 

Independent of the results of sections 2, 3 and 5, the important new information 
concerning the relationship between the particle's momentum and velocity according to 
the MPD equations at condition (4) is presented in expressions (37)-(40). 

In another paper we plan to carry out a more detailed analysis of the possible 
choice of values E, J and to present the results of a complex investigation of the highly 
relativistic motions of a spinning particle in the Kerr space-time according to the exact 
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MPD equations (22), (A.3)-(A.9). 



Appendix. Explicit form of the seven equations described in section 3 

First, we stress that the complex exphcit form of the seven equations in a Kerr 
metric, which were discussed in section 3, is determined by the much more complicated 
expressions for the components of the Riemann tensor and the Christoffel symbols for 
this metric than for Schwarzschild's one. We do not write these expressions here for 
brevity because they are presented in other papers, for example, in the Appendix of [8] . 
We need introduce the notation (in addition to (18)-(21)): 

z ^ yl + cos^ 1/2, q = yi{yi - 2) + , ip = yl - cos^ y2, 

Tj^Syf- cos^ ^2, x = 2/i + i^y\- 'io? cos^ ^2, (^-1) 

where a is equal to ajM. One can see that very quantities analogous to (A.l) are 
presented in the expressions for many components of the Riemann tensor and the 
Christoffel symbols [8]. 

To achieve more compact form of the equations we use other notation as well: 

Pi = -zy^q~^, P2 = -zye, 

P3 = [2ayiy8 sin^ y2 - yvizx + ^a^Vi sin^ ^2) sin^ Z/2]^"\ 
P4 = [2ayxy'jS\T?y2 + ys{z - 2yx)\z~'^; 
ci = yiyw sin y2-{,z- 2yi)yQyiiq'^ sin"^ ?/2, 
C2 = ybVii^z - 2yi)g-^sin-^y2 - Z/7Z/9?siny2, 

C3 = {y&y^q - ysyio) siny2; 

di = -{2aq~^yiyey 11 + ysVio) siny2, 
d2 = {2aq'^yiy^yii + gt/sl/g) sin 7/2, 

4 = {ybVio - qy&y9) sin 1/2; 



p = 2ayiy'j sin^ y2 + {z- 2yi)ys. {A.2) 
Then the pointed out in section 3 all seven equations can be written as: 

y^ijb + Z/ioye + yum = A- y^Qi - yioQ2 - yiiQs, {A.3) 

Piijb + ^2^6 + P2.y7 + Pirn = 

= -PiQi - P2Q2 - P3Q3 - PiQi, (AA) 

ciyb + C2ye + c^m = C - ciQi - C2Q2 - C3Q3 + E, (A. 5) 
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diy5 + d2y6 + dsys = -D - diQi - d2Q2 - dsQ^ - J, {A.6) 

py9 = oty^y7y%{.a^ - yl)q~^ sin^ 2/2 - a^yiy&y-iy^z'^ sin^ 7/2 sin 27/2 + a'^y^ysysiyi - 1)?"^ 

X sin^ y2 - O.ba^yeysygiz - 2yi)z''^ sin2?/2 + ayfyni/jq'^ + 2ayiy^yn{yiz{z - 2yi) 
-a^i) sin^ 7/2)^"^?"^ sin yl + aylyu {z - 2yi)tpz'^q~^ + ayiyr^y^yioq'^ sin 2y2 + O.ba'^y^ 
xysyioiz - Ayi)z~^q~^ sin 2y2 - 2ayiy^yQyiiq~^ coty2 + 2aylyeyryioz'^ sin^ 7/2 + ymys 
xyio{z - 2yi)z'^ + yvysyuiyiziz - 2yif - a'^'ip{z - 4yi) sin^ y2)z~'^q~^ + j3[2ayly7{q{z 
-3yf) + yiz{yi - sin^ 7/2 + ylysq''^ {^qi^ + a'^z{l - yi) sin^ 7/2) + 2ayiy5y6y7 

x{z + 2a^ sin^ y2)q'^ sin 2?/2 + a^y^yeysiz - 4yi)g"^ sin 2|/2 + 2aylyly7 sin^ 1/2 
+yiyly?,{.z - 2yi) - yl{z - 2yi)^ljz~^ + 2ayjylilj{z - 3yi)z~'^ sin^ 7/2 + 2ayiy^[yizxiz 
-2yi) + sin^ y2{zq - 2o?y\ sin^ 7/2 + ^y\)\z~'^q~^ sin^ t/2 + 2ayxz{yjy^ - y->y-j)q~^ sin^ 7/2 
^z{z - 2yi){y5ys - y^ybjO."^ + yyy%{yiz^{z - 2yi) - a^ilj{z - 6yi) sin^ y2)z~'^ sin^ ^2], (^-7) 

pyio = ay5y6yn{2yf - z)q'^ - 2a^yiyr,yjy^z~^ sin^ 7/2 cos 7/2 + O.Sa^ysysyglSyi - 2^)2;"^ 

X sin 2^/2 + 2ayly5y7yioZ~^ sin^ I/2 + l/il/5l/8l/io(^ - 2i/i)^"^ - 2ayiylynCoty2 
+ayiyjyu{z + 2a^yi sin^ 2/2)^"^ sin 2^2 + 2ayiylyii{2yi - z)z~'^coty2 + ay^y-jy^q 
x{z- ^yl)z~^ sin^ 7/2 + y&ysy^qi^yl - z{yi + l))^;"^ + ayiy^y-ryw sin 2^2 - O.Sa^yeysyio 
X sin 2y2 + yrysyiiiz^ + 2yiz{a^ sin^ 7/2 - ^) - 8q;^?/i sin^ Z/2)^~^ cot 7/2 + P[^.^a^ylys{2yi 
-z)q~^ sin 2^2 - 2a^yiylyjq-^ sin^ 1/2 cos t/2 - 2ay5yey77] sin^ 1/2 + 2y5yQy8{Ayl - z{yi + 1)) 
+ayiylyj{2z + 3q;^ sin^ 7/2) sin2y2 + a^ylys{z - Qyi) sint/2 cos 7/2 + 2ayiz{y6y7 - yTy^) 
X sin^ y2 + z{z - 2yi){yey8 - ysye) + ot^y\yl{z - 2yi)z~'^ sin 2^2 + 2ayiy^{xz'^ 
+Aa'^yiz sin^ 7/2 + 2a\ sin^ Z/2)^~^ sin^ 7/2 cos 7/2 + 2ayiy-jyl{a^yi svc? 7/2 - - 2yi))z'''^ 

X sin 2|/2 + 0.5y^ys{z^q + 2a^yi{zq - 6yix) sin^ 2/2)^"^ sin 27/2], (A. 8) 

pyii = -2a^yiy6y7yiiz~^ sin^ 7/2 cos 7/2 + ay^yiyii [{^yl - z)x - 2a^yl sit? 7/2 - 

X sin^ y2 + ysysZ/ii [2/1(2 - 2yi)2 + a^{z- 2yf) sin^ 2/2]^"^?"^ + yeysyuiziz - 2yi) + 2a^yi 
X sin^ y2\z~'^ cot ^2 - ay^ygqia^i/j + yf{z + 2yl))z~^ sin^ + y7y8y9q{zyi{2yi - z) + 
+q;^ sin^ 2/2))^"^ sin^ 2/2 - aylygq^pz'^ sin^ 7/2 + 2a^yiyjyioqz~'^ sin^ 2/2 cos 2/2 - 0.52/72/82/10? 
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X (2;^ + Aa'^yi sin^ y2)z~'^ sin 2^2 + ayiylyioqz'^ sin 2^2 + /3[-2ay5yy{xip + 2yfz)z~^ sin^ ^2 
+4:a^yiyfiyjqz~^ sin^ ^2 cos ?/2 - 2ay^yl'il)z~'^ sin^ ^2 + 2ayiyeylqz~^ sin 2^/2 + ^g(z/7y8 
-I/sJ^t) sin^ 1/2 - y&yiy&zq sin 2|/2 - 2yr,yjys{yiz{z - 2yi) - ijjix + c? sin^ 2/2))^"^ sin^ I/2], (A.9) 
where 

^ = -2ay5y6y9yior]z~^ cos 1/2 - 2ay^yjy<;^yiiriz'^{x + sin^ 2/2) cos 2/2 
-Qayiyayrygyuq^z-^ sin ^2 + ^ymy%y%yi\qiz~'^ sin"^ ^2 + ^Oiyiy^ygyioq^x^-^ sin^ ^2 
-^ymysygym^z'^ix + sin^ ^2) sin y2 + Qoiyiylygyioq^z"^ sin ?/2 
+2ai/6l/7yioyii^2;~'^(2x + sin^ t/a) cos?/2 - Qayiy^yiyimiiXQ.'^ sin?/2 
+6a\i|/5y8yioyiiCg"^2;"'' sin 1/2 + ayly^qrjz'^ cosj/2 
+Oiy^ygr)q{x^ + 2qa^ sin^ ^2)2;"^ sin^ ^2 cos y2 
-2o?y'jy^ylr]q{Zx - ^yi)z~^ sin^ ^2 cos + ar??(2g + sin^ y2)z~^ylyl cos ^2 

+avyhw'l~^^~^ cos 1/2 - tt^?(2x^ + tt^g sin^ l/2)l/?l/io^"^ sin^ I/2 cos 1/2 
-Q;r/(g + 2a2 sin^ y2)ylylQZ~^ cos 1/2 + 2a^?7(3x - 2yi) x yiy^yloZ'^ sin^ I/2 cos 7/2 
+Q;77(g + 2a^ sin^ y2)ylyliq~^ sin"^ ^2 cos ^2 - Q;77(2g + sin^ ^2) 
^ylyliff^z'^ sin~^ y2 cos ^2 + Goi^yiysyeyuq'^ z"^ sin~^ y2 
-4.a^ylr]yjyl^q~^z~^ s\v? y2 cos ^2 
-2a^yir]{q - sin^ I/2) (1 + sin^ ?/2)|/7|/8?/iig"^^~^ cos 7/2 

-q;77(^ - 2|/i) (g - q;^ sin^ y2)ylyliq~^ z~^ cos ^2; (^.10) 

Qi = (yi? - ^(yi - l))y5^"^?"^ - qyiylz'^ - qiyiz^ - a'^ipsin^y2)yjz~^sm'^y2 
+qi/jylz~^ - a^yr,yQZ~'^ sin 2^2 - 2aqijjy7y8 z~^ sin^ ^2, 
= -O.Sa^gZ"^ sin2?/2 + 0.5a'^ylz~'^q~'^ sin2|/2 - 0.5yj{z^x 
+2Q!^|/i(x + ^) sin^ 1/2)^"^ sin 2y2 - a^yiylz'^ sin 2y2 + 2yiy^yez~^ + 2ayiy7ysxz~^ sin 2^/2, 
<53 = 2y^y'j{yiz{z - 2yi) - a'^^ sin^ y 2) z~'^q~^ + 2ay5y8ipz~'^q~^ 
+2yey-r{z'^ + 2a^yi sin^ y2)z~'^ cot y2 - Aay^y^y^z''^ cot ys- 
(54 = -2ay5y7{2yfz + ^x)^:"^?"^ sin^ I/2 + 2yr,y8il)xz~'^ q~^ 

+2a^yiyQy7Z~'^ sin^ 7/2 sin 2^/2 - 2a^yiyQysZ~^ sin 2^/2; (^-H) 

C = -(1 - 2yiz~^)y8 - 2ayiy7Z~^ sin^ 2/2 
+2a'^qyiy7ygz''^ sin^ y2 cos - 2aqyiy8ygz''^ cos ^2 
+xV'y7yio2;~^ sin y2 - aysyioV'^"^ sin y2 
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-2a'^yiy5yiiq~^z~'^ cos y2 - yeVu'ipz''^ sin"^ (^-12) 

D = -2a^yxy-jy^z~^ sin^ ^2 cos y2 + qiz^ + 2Q;^yi sin^ ^2) 
xysVdZ'^ cos 1/2 - ayjy loixi^ + 2ylz)z~^ s\v? y2 

-ymo{yiz{z - 2yi) - a^ipsw?y2)z~^siny2 
+2ayiy5yuxq~^z~'^ cosy2 + aipyeyuz''^ siny2 
+p^^[2ayiyj{zx + 2a;^yi sin^ Z/2)^~^ sin^ ^2 + qzyjVs 
x{l - Sa'^yfq'^z-'' sm'' y2) - 2ayi{z - 2y^)z-'yl] sin'' y2; (A 13) 



/3 = 1/51/9 + 1/61/10 + 1/71/11- (^-14) 

We stress that both equations (A.3)-(A.9) and expressions (A.10)-(A.13) become 
much simpler in the case of the Schwarzschild space-time, when a = a/M = 0. Then, 
for example, instead of long expression (A. 10) we have 

A = dyiyeysygyuq^z'^ sin~^ y2 - ^yiViVsy^ymiz'^ {x + sin^ ^2) siny2. (^.15) 
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